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Can a primordial black hole or wormhole grow as fast as the universe? 
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This review addresses the issue of whether there are physically realistic self-similar solutions in 
which a primordial black hole is attached to an exact or asymptotically Friedmann model for an 
equation of state of the form ^ = (7 — l)pc^. In the positive pressure case (1 < 7 < 2), there is no 
such solution when the black hole is attached to an exact Friedmann background via a sonic point. 
However, it has been claimed that there is a one-parameter family of asymptotically Friedmann 
black hole solutions providing the ratio of the black hole size to the cosmological horizon size is 
in a narrow range above some critical value. There are also "universal" black holes in which the 
black hole has an apparent horizon but no event horizon. It turns out that both these types of 
solution are only asymptotically gi^asz-Friedmann, because they contain a solid angle deficit at large 
distances, but they are not necessarily excluded observationally. Such solutions may also exist in 
the 2/3<7<<l case, although this has not been demonstrated explicitly. In the stiff case (7 = 2), 
there is no self-similar solution in an exact background unless the matter turns into null dust before 
entering the event horizon, which is a contrived and probably unphysical situation. However, there 
may be asymptotically quasi-Friedmann solutions without a sonic point which contain universal 
black holes. In the negative pressure case (0 < 7 < 2/3), corresponding to a dark-energy-dominated 
universe, there is a one-parameter family of black hole solutions which are properly asymptotically 
Friedmann (in the sense that there is no angle deficit) and such solutions may arise naturally in 
the inflationary scenario. The ratio of the black hole size to the cosmological horizon size must now 
be below some critical value, so the range is more extended than in the positive pressure case and 
one needs less fine-tuning. If one tries to make a black hole which is larger than this, one finds a 
self-similar solution which connects two asymptotic regions, one being exactly Friedmann and the 
other asymptotically quasi-Friedmann. This might be regarded as a cosmological wormhole solution 
providing one defines a wormhole throat quasi-locally in terms of a non-vanishing minimal area 
on a spacelike hypersurface. The possibility of self-similar black holes in phantom fluids (7 < 0), 
where the black hole shrinks as the big rip singularity is approached, or tachyonic fluids (7 > 2) 
remains unclear. We also consider the possibility of self-similar black hole solutions in a universe 
dominated by a scalar fleld. If the fleld is massless, the situation resembles the stiff fluid case, so 
any black hole solution is again contrived, although there may still be universal black hole solutions. 
The situation is less clear if the scalar field is rolling down a potential and therefore massive, as 
in the quintessence scenario. Although no explicit asymptotically Friedmann black hole solutions 
of this kind are known, they are not excluded and comparison with the < 7 < 2/3 perfect fluid 
case suggests that they should exist if the black hole is not too large. This implies that a black 
hole might grow as fast as the cosmological horizon in a quintessence-dominated universe in some 
circumstances, supporting the proposal that accretion onto primordial black holes may have played 
a role in the production of the supermassive black holes in galactic nuclei. 

PACS numbers: 04.70.Bw, 97.60.Lf, 04.40.Nr, 04.25.Dm, 95.35.-hd 



I. HISTORICAL INTRODUCTION 

Over the last 40 years there has been much interest 
in how fast a black hole formed in the early universe, 
when the density is usually radiation-dominated, would 
grow. As first pointed out by Zel'dovich and Novikov [1], 
a simple Bondi-type accretion analysis suggests that a 
primordial black hole (PBH) would not grow much at all 
if it were much smaller than the cosmological horizon at 
formation but that it could grow at the same rate as the 
universe if its initial size were comparable to it. (The 
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term "cosmological horizon" should here be interpreted 
as the Hubble horizon if the PBHs form after an infla- 
tionary period but the particle horizon otherwise.) One 
might expect the latter situation to apply, since a PBH 
must be bigger than the Jeans length at formation [2], 
so this suggests that any PBH might grow to the horizon 
mass at the end of the radiation era, which is around 
10^^ Mq. Since there is no evidence for such enormous 
black holes, for a while it was assumed that no PBHs 
ever formed. 

However, the validity of the Zel'dovich- Novikov calcu- 
lation is questionable when the black hole size is compa- 
rable to the horizon size because it neglects the expansion 
of the universe and is not fully relativistic. Indeed, the 
conclusion that a PBH could grow at the same rate as the 
universe in the radiation-dominated era was disproved 
by Carr and Hawking [3]. They demonstrated this by 
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proving that there is no self-similar solution which con- 
tains a black hole attached to an exact flat Friedmann 
background via a sonic point (i.e. in which the black 
hole forms by purely causal processes). The Zel'dovich- 
Novikov prediction is therefore definitely misleading in 
this case. Since the PBH must soon fall well below the 
horizon size, when their argument should be valid, this 
suggests that PBHs would not grow much at all. 

This gave the subject of PBHs a new lease of life and 
motivated Hawking to consider the quantum effects asso- 
ciated with black holes (since only PBHs could be small 
enough for these to be significant). Ultimately, this led 
to his discovery of black hole radiation [4] , so it is ironic 
that a consideration of PBH accretion led to the conclu- 
sion that they evaporate! 

Carr and Hawking also claimed that there are self- 
similar solutions which are asymptotically - rather than 
exactly - Friedmann at large distances from the black 
hole. However, these correspond to special acausal ini- 
tial conditions, in which matter is effectively thrown into 
the black hole at every distance; they do not contain a 
sonic point because they are supersonic everywhere. In- 
deed, such solutions exist in the "dust" case, when the 
cosmological fluid is pressureless [3]. They even found 
solutions in which the whole universe is in some sense in- 
side a black hole; these are now termed "universal" black 
holes. 

Subsequently, the Carr-Hawking analysis was extended 
to perfect fluids with equation of state p = — ^)pc^ 
where p is the pressure, p is the mass density and 7 is 
a constant (4/3 in the radiation case). This is the most 
general form for a barotropic equation of state compatible 
with self-similarity. It was proved that the non-existence 
of self- similar black holes in an exact Friedmann back- 
ground applies for all values of 7 in the range 1 to 2 [5, 6]. 
Indeed, the only physical self-similar solution which can 
be attached to an exact external Friedmann solution via 
a sonic point is Friedmann itself; as the radial coordinate 
decreases, the other solutions either enter a negative mass 
regime or encounter another sonic point where the pres- 
sure gradient diverges [6]. As in the radiation case, there 
are still acausal black hole solutions but these are again 
supersonic everywhere; the asymptotically Friedmann so- 
lutions which reach a sonic point represent density per- 
turbations which grow at the same rate as the universe 
rather than black holes [7]. 

Later it was realized that none of these positive- 
pressure self-similar solutions are strictly asymptotically 
Friedmann after all [8]: there is a solid angle deficit 
at large distances which might in principle show up 
in the angular diameter test. It would therefore be 
more accurate to describe them as asymptotically "quasi- 
Friedmann". Such solutions are not excluded observa- 
tionally, at least for some parameter range, but they are 
not physically well- motivated. 

The attempt to extend the analysis to stiff fluids 
(7 = 2) led to some controversy. Lin et al. [9] at first 
claimed that there is a self-similar black hole solution in 
an exact Friedmann background in this case. However, 
Bicknell and Henriksen [10] showed that this conclusion 



is invalid because Lin et al. had misidentified the point 
corresponding to the black hole event horizon. Bicknell 
and Henriksen did manage to construct a numerical self- 
similar solution containing a black hole but it required 
the stiff fluid to turn into null dust at some point. Al- 
though this might seem rather contrived. Reed and Hen- 
riksen [11] later found a solution of this kind by general- 
izing some work of Hacyan [12], involving a self-similar 
Vaidya model. However, even this model now seems im- 
plausible. The only possibility might be universal black 
hole solutions which are asymptotically quasi- Friedmann. 

It is also interesting to consider the growth of a black 
hole when the density of the universe is dominated by a 
scalar field, as expected in many cosmological contexts. 
If the scalar field is massless (i.e. if there is no scalar po- 
tential), then one might expect the same result to apply 
as in the stiff fluid analysis, since it is well known that a 
scalar field is equivalent to a stiff fluid provided the gra- 
dient is everywhere timelike [13]. Indeed, the conclusion 
that there is no self-similar non-universal black hole solu- 
tion in an exact or asymptotically Friedmann background 
dominated by a massless scalar field is supported by both 
numerical studies [14] and analytical calculations [15]. 

However, the situation is more complicated if there is 
a scalar potential (i.e. if the scalar field is massive, the 
mass being associated with second derivative of the po- 
tential) and the discovery that the universe is currently 
accelerating suggests that this may be the case at the 
present epoch [17]. This has led to a study of black hole 
accretion in quintessence- dominated universes. Indeed, 
an argument similar to that advocated by Zel'dovich and 
Novikov has resurfaced in this context in order to ex- 
plain the origin of the 10^ to 10^ Mq black holes thought 
to reside in galactic nuclei [18]. While there are sev- 
eral scenarios for the formation of such supermassive 
black holes, one possibility is that they originated in 
the early universe and grew self-similarly through ac- 
cretion of quintessence before cosmological nucleosynthe- 
sis [19, 20], as well as by purely astrophysical processes 
later. 

Since this proposal is motivated by a Bondi-type argu- 
ment which neglects the cosmological expansion, it is just 
as questionable as the original Zel'dovich-Novikov one. 
This has led to the search for self-similar black hole solu- 
tions in quintessence-dominated universes. In this case, 
the similarity assumption requires that the potential have 
an exponential form. Our analysis in Ref. [15] then shows 
that there is no self-similar solution with a black hole in 
an exact or asymptotically Friedmann or asymptotically 
quasi-Friedmann background if the universe is decelerat- 
ing and no such solution in an exact Friedmann back- 
ground if it is accelerating. However, this does not prove 
non-existence in an asymtotically Friedmann or quasi- 
Friedmann background for the case in which the back- 
ground is accelerating. Kyo et al. [21] have shown that 
there is a one-parameter family of self-similar asymptot- 
ically Friedmann solutions in this case, although it is un- 
clear whether they can contain black holes. 

The acceleration of the universe can also be explained 
if its density is dominated by a perfect fluid with < 7 < 
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2/3, so this has motivated us to look for self-similar black 
hole solutions in this case [22, 23]. Such fluids are very 
different from positive-pressure ones, since there are no 
sound-waves (the sound-speed \fpjp being imaginary), 
so one might expect the conclusion about self-similar 
black hole solutions to be different. We describe such a 
fluid as "dark energy", although this term is sometimes 
used more generally and may indeed include quintessence 
itself. One might regard quintessence as a form of dark 
energy in which the parameter 7, rather than being con- 
stant, may vary as the scalar field rolls down its potential. 
However, there is an important physical difference be- 
tween (constant 7) dark energy and quintessence because 
there are sound-waves in the latter case, the sound-speed 
being the speed of light at short wavelengths [15]. 

Since there are no sound-waves for a dark energy 
fluid, there can be no black hole solutions in an ex- 
act Friedmann background. However, as discussed in 
Refs. [22, 23], there are asymptotically Friedmann solu- 
tions containing black holes in this case. Indeed, unlike 
the positive-pressure case, these are genuinely asymp- 
totically Friedmann rather than asymptotically quasi- 
Friedmann and the associated inhomogeneities may arise 
naturally in the inflationary scenario. There is a one- 
parameter family of such solutions but (in contrast to 
the implication of the Zel'dovich-Novikov argument) they 
only exist if the black hole is not too large compared to 
the particle horizon. These solutions are not analogous 
to the positive-pressure solutions which Carr and Hawk- 
ing were originally seeking but they might nevertheless 
be physical. 

If one tries to find an asymptotically Friedmann self- 
similar solution with a black hole which is larger than 
the upper limit in the < 7 < 2/3 situation, one obtains 
a cosmological wormhole instead. The transition occurs 
as the black hole apparent horizon approaches the cos- 
mological apparent horizon, after which both horizons 
disappear. This is in contrast to the 2/3 < 7 < 2 case, 
where the two apparent horizons never merge [5] and one 
tends to a separate closed universe as the black hole size 
increases. (However, the separate universe case is not 
itself self-similar.) In the wormhole solution the metric 
tends to an asymptotically K ant owski- Sachs form as one 
approaches the wormhole throat, this corresponding to a 
minimum physical radius, and the solution then connects 
to another asymptotically Friedmannn or asymptotically 
quasi-Friedmann universe. 

This paper provides a comprehensive discussion of all 
these solutions. For the most part, we will avoid math- 
ematical technicalities, so the number of equations is 
minimized. Although the paper is intended as a gen- 
eral review of previous work in this area, we believe 
that bringing all the cases together is illuminating and 
leads to some original insights. In successive sections, we 
consider positive-pressure fluids, stiff fluids, scalar and 
quintessence fields, dark-energy fluids and finally more 
exotic possibilities (phantom fluids with 7 < 0, negative 
pressure fluids with 2/3 < 7 < 1 and tachyonic fluids 
with 7 > 2). We conclude that there are certainly self- 
similar asymptotically Friedmann black hole solutions in 



the < 7 < 2/3 case and there may also be in the 
quintessence case, so it is interesting that these situations 
may be observationally favoured at the present epoch. In 
the final section we draw some general conclusions. Two 
appendices clarify the dimensions of various quantities 
used in our analysis and the connection between relevant 
energy conditions. 

It should be stressed that this paper is not intended to 
be a review of the more general problem of black holes 
in a cosmological background, although that problem is 
of great interest in its own right [24-29]. It is also much 
more narrowly focussed than the earlier review of self- 
similar solutions in general relativity by Carr and Co- 
ley [30]. 



II. POSITIVE-PRESSURE FLUID CASE 

The Zel'dovich-Novikov argument [1] is based on a 
Bondi [16] analysis, in which one neglects the cosmo- 
logical expansion and assumes that the background fluid 
will be swallowed from within an accretion radius Ra = 
2GMj(?g^ where Cg is the sound-speed and M is the back 
hole mass. The accretion rate for the black hole is then 
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where p is the mass density at the accretion radius, which 
is taken to be the background cosmological density. In 
the radiation-dominated era, which was the case consid- 
ered by Zel'dovich and Novikov, Cg = c/\/3 and so Ra 
is just three times the Schwarzschild radius of the black 
hole. Since the mass density in a Friedmann radiation- 
dominated model is p = 3/(327rGt^), where t is the cos- 
mological time from the big bang, Eq. (2.1) then gives 
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Note that this equation cannot be precisely correct. As 
discussed later, the density near the black hole is not ex- 
actly the background value and the accretion radius is 
not exactly 2GM/c^, but these effects merely modify the 
effective value of K. Equation (2.2) can now be inte- 
grated to give 
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where Mq is the mass of the PBH at its formation time to 
and K~^to is roughly the mass within the cosmological 
horizon at that time. (Both the particle horizon and 
Hubble horizon mass are exactly c^t/G in the radiation 
era, which is a factor 9\/3/2 ^ 8 larger.) For Mq <^ 
K~^to, Eq. (2.3) implies negligible accretion. However, 
for Mq — K~^tQ, it predicts M = K~^t, in which case the 
black hole always has a size comparable to the particle 
horizon. For Mq > K~^tQ, the mass diverges at a time 



to 1 



to 



KMn 



(2.4) 



4 



(i.e. at a time comparable to to unless Mq is very close 
to K~^to). These three behaviours are illustrated by the 
upper frame in Fig. 1. 
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FIG. 1: Schematic figure showing the growth of PBHs with 
positive pressure for a Bondi-type analysis which neglects the 
cosmic expansion (top) and a full relativistic analysis which 
allows for it (bottom). 



The middle case corresponds to self-similar growth and 
is precisely the situation which one might expect to per- 
tain, since an overdense region needs to be larger than 
the Jeans length at formation but not so large as to form 
a separate closed universe [31], both these length scales 
being of order the cosmological horizon size. This sug- 
gests that a PBH might in principle continue to grow 
as fast as the cosmological horizon until the end of the 
radiation-dominated era, when its mass would be of order 
IO^^Mq. The existence of many such holes is excluded, 
for example, by observations of the large-sale bulk flows 
of galaxies [32] , which might lead one to doubt that PBHs 
ever formed. Of course, the self-similar solution requires 
very fine-tuning of the black hole mass in the Zel'dovich- 
Novikov argument but, providing the black holes span a 
spectrum of masses, one would expect at least some of 
them to satisfy this condition. 

That the behaviours shown the upper frame in Fig. 1 
cannot be correct in a more precise calculation is indi- 
cated by the lower frame in Fig. 1. This shows the re- 
sults of fully relativistic numerical calculations for the 
same initial black hole masses as in the upper frame. We 
next investigate why the simple Bondi analysis fails and 
how it can be improved. 



A. Improved analysis neglecting cosmic expansion 

The Zel'dovich-Novikov argument can be refined in 
various ways. In particular, Harada and Carr [34] im- 
proved the analysis by allowing the fluid to have a more 
general barotropic equation of state p = (7 — l)pc^ with 
1 < 7 < 2 and by including a relativistic focussing fac- 
tor a in Eq. (2.1). In this case, the cosmological density 
is p = l/(67rG7^t^), so Eq. (2.3) remains valid but the 
factor K becomes 



K 



SGa 



3c3(7- 1)3/2^2 



(2.5) 



(Their expression was actually a factor of 4 smaller than 
this because they took the accretion radius to be GM/c^ 
rather than 2GM/(?^.) The length-scale associated with 
the self-similar black hole solution is 
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(or a factor 4 larger in the Harada-Carr analysis), where 
H = 2/(37t) is the Hubble parameter. They also showed 
that an overdense region is a separate closed universe 
rather than part of our universe if its size exceeds 
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This expression applies for 7 > 2/3 and not just for 
1 < 7 < 2. The self-similar scale is always less than 
the separate-unverse scale in the present analysis but it 
goes above the self-similar scale for 7 > 1.6 in the orig- 
inal Harada-Carr analysis. They inferred that - even 
before attempting a more exact calculation - one should 
have reservations about the self-similar prediction for suf- 
ficiently hard equations of state. 

Harada and Carr gave precise expressions for all rele- 
vant length scales in terms of the equation of state pa- 
rameter and the functional forms are indicated in Fig. 2. 
The other scales shown there are the Jeans length [33] 
and the particle horizon size: 
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Since there is some uncertainty in the appropriate accre- 
tion radius to use, various possibilities for the self-similar 
scale are shown; these all take the focussing factor a to 
1. Note that some of the scales in Fig. 2 can be extended 
down to 7 = 2/3 and they can all be extended above 

7 = 2. 

Babichev et al. [35] have further refined the accretion 
analysis by replacing Eq. (2.1) with 
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where the subscript oo indicates asymptotic values and 
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and this is also shown in Fig. 2. Note that ah three 
expression for Rs shown there scale as (7 — !)^/^ as 7 ^ 1 
and the Babichev et al. expression actually coincides 
with the original Harada-Carr one in this limit. The ratio 
HRs/c is 0.5 at 7 = 2 and asymptotes to 8/(3>/3) ^1.6 
as 7 ^ 00, so the self-similar scale is always smaller than 
the separate-universe scale. 

Clearly these equations do not apply for 7 < 1 since 
the value of A given by Eq. (2.10) is not well defined 
then. However, Babichev et al. argue that one can take 
A = 4 in this situation, corresponding to the stiff fiuid 
value. In this case, Eq. (2.9) gives a sensible accretion 
rate, even though the accretion radius given by Eq. (2.12) 
is not well defined. We will use this result in the later 
discussion. 



FIG. 2: The length scales associated with the Jeans condi- 
tion, the particle horizon, the Hubble horizon, the separate- 
universe condition and the self-similar black hole solution are 
plotted against /c = 7 — 1. Three self-similar scales are shown, 
corresponding to the accretion radius being (1) GM/cg, (2) 
2GM/cl and (3) the expression given by Eq. (2.12). 



This allows for the fact that the density of the fiuid at 
the accretion radius is not the asymptotic cosmological 
value and it also incorporates the relativistic correction 
to Eq. (2.2) associated with the pressure. The Babichev 
et al. analysis is not confined to the p = (7 — l)pc^ 
equation of state, but we assume this here because of 
the self-similarity requirement. Thus we again obtain 
Eq. (2.2) but with K is replaced by 
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SO the mass for self-similar growth now has a more com- 
plicated dependence on the equation of state parameter 
than in the Harada-Carr analysis. As 7 decreases from 2 
to 1, Eq. (2.10) indicates that the factor A increases from 
4 to 00, passing through 6a/3 at 7 = 4/3 and diverging 
like [e/(7-l)]^/^/4 as 7 ^ 1. Comparison with Eq. (2.1) 
also shows that the "effective" accretion radius is 
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The last expression shows the deviation from the naive 
expression Ra = 2GM/cl. For example, this is a factor 
\/2 for a stiff fluid and a/2/3 for a radiation fluid. Ra 
diverges in the limit 7^-1 and has the form 
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the numerical coefficient being almost exactly 1. The 
length-scale associated with the self-similar solution is 
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B. Improved analysis allowing for cosmic expansion 

In all the above analyses, the existence of the self- 
similar black hole solution relates to the fact that the 
density scales as p ex in a ffat Friedmann background. 
However, none of the refinements discussed above allows 
for the cosmological expansion itself, which one might 
expect to oppose accretion, so one should be wary of the 
self-similar prediction. 

Carr and Hawking's refutation of the Zel'dovich- 
Novikov argument was based on the study of spheri- 
cally symmetric self-similar solutions to Einstein's equa- 
tions [3]. Such solutions have the characteristic that all 
dimensionless quantities depend only on the "similarity" 
variable z = r/t, where r is the comoving radial coor- 
dinate and t is the cosmological time [36]. The metric 
takes the form 

ds^ = -e^^^'Ut^ + e'^'^'Ur^ + r^S{zfd^\ (2.15) 

where d^t^ = dO'^ + sin^ Od^p^ and we now put c = 1 
throughout this paper unless specified otherwise. A key 
role is played by the "velocity" function 
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which represents the fiuid velocity relative to the simi- 
larity surfaces of constant z (which are expanding or col- 
lapsing spheres). Providing a certain regularity condition 
is satisfied [15], a value of z where V = 1 corresponds to 
a null surface and is termed a "similarity horizon" . The 
fiat Friedmann solution is itself self-similar and has 
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The point where V = 1 corresponds to the cosmologi- 
cal particle horizon for 7 > 2/3. (It corresponds to a 
cosmological event horizon for 7 < 2/3 but we assume 
1 < 7 < 2 in this section.) If one had a black hole in 
a Friedmann background, there would be two values of 
z for which V = 1, the inner one being associated with 
the black hole event horizon and the outer one with the 
cosmological particle horizon. 

Values of z for which V = ^/dipjdp = y^7 — 1 are 
also important because this may correspond to a sonic 
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point (i.e. the pressure gradient may be discontinuous). 
For the radiation-dominated universe considered by Carr 
and Hawking, the sound speed is F = 1/a/3. Since the 
Friedmann solution itself has such a point, one might 
envisage a self-similar solution in which one attaches a 
black hole interior to an exact Friedmann exterior via 
a sound-wave. If such a solution existed, it would repre- 
sent a black hole growing at the same rate as the universe 
and V{z) would then have the form indicated by the bro- 
ken curve in Fig. 3. However, Carr and Hawking showed 
that such a solution is not possible [3] . This implies that 
a black hole formed by local processes can never grow 
as fast as the universe, even if it starts off with a size 
comparable to the particle horizon. Thus the Zel'dovich- 
Novikov prediction cannot be correct in this case. 

V 



universal black hole 




FIG. 3: Schematic figure of the function V(z) for self-similar 
asymptotically quasi-Friedmann black hole solutions, showing 
the two similarity horizons. The solutions are described by 
a single parameter and the minimum of V reaches the sonic 
value for solution. A black hole connected to an exact Fried- 
mann background via a sound-wave would look like the dotted 
curve but such solutions cannot exist. 

This conclusion was soon supported by detailed hydro- 
dynamical calculations. If the formation and evolution 
of a PBH in a positive-pressure fluid with a local per- 
turbation is simulated numerically, it is found that the 
PBH soon becomes much smaller than the cosmological 
horizon [37-39]. In this situation, one would expect the 
Zel'dovich-Novikov argument to apply, so one does not 
anticipate the black hole growing very much at all. This 
also applies for holes which are initially larger than the 
horizon [14], even though their mass should quickly di- 
verge according to the Bondi-type analysis. 

On the other hand, Carr and Hawking also constructed 
a one-parameter family of spherically symmetric self- 
similar PBH solutions for a radiation (7 = 4/3) or dust 
(7 = 1) fluid which they described as asymptotically 
Friedmann [3]. In these solutions, as z decreases, the 
velocity V decreases to a minimum between 1 and the 
sound-speed l/VS. V then rises through 1 at the black 
hole event horizon to infinity at the black hole singular- 
ity (corresponding to a finite value of z but zero physical 
distance). Two such solutions are shown in Fig. 3; we 



describe these as "acausal" since they are supersonic ev- 
erywhere and the perturbation extends beyond the cos- 
mological horizon. There are also solutions in which the 
minimum of V exceeds 1, two examples of these being in- 
cluded. This corresponds to a "universal" black hole [40] 
in which there is no black hole or cosmological particle 
horizon. In this case, the entire universe might be re- 
garded as being inside the black hole, although there are 
still black hole and cosmological apparent horizons [5]. 
The conclusion that there are self-similar black hole so- 
lutions in a universe which is asymptotically Friedmann 
but no such solutions in an exactly Friedmann universe 
was soon generalized to perfect fluids with equation of 
state p = (7 - l)p with 1 < 7 < 2 [5, 6]. 

It is interesting to examine these asymptotically Fried- 
mann solutions in more detail. The single parameter Bq 
which describes them is related to the perturbation at 
large z of the factor S{z) which appears in metric (2.15): 

S = SF{z)e''^'\ B^Bo^B^z-^^^^-^y^^^\ (2.18) 

where Bi is determined by Bq and the value of 7. This 
is also related to the density perturbation at large dis- 
tances: 

W = SnGpr^ = W^{z)e^^'\ A ^ Aiz-^^^^-^^^^^^\ 

(2.19) 

where Ai = {2 — 37)^1. For a sufficiently negative value 
of 5 corresponding to a sufficiently large overdensity, 
V will reach a minimum as z decreases. This minimum 
will decrease as Bq increases but it will exceed 1 for Bq < 
/3i < 0, where the value /3i depends on 7. These solutions 
correspond to universal black holes. 

On the other hand, the minimum will lie between V = 
1 and for Pi < Bq < < 0, corresponding to 

the black hole solutions discovered by Carr and Hawking. 
Here ^2 also depends on 7 and is the value of Bq for which 
the minimum of V reaches y^7 — 1 [7] . Both these types 
of solution are also illustrated in Fig. 3. Note that the 
ratio of the black hole event horizon size to the particle 
horizon size decreases as Bq increases and so must always 
exceed the value associated with the solution with Fmin = 
— 1. It might seem plausible to identify this with the 
Jeans criterion for gravitational collapse, in which case 
the ratio of the black event horizon and particle horizon 
should also have a minimum of order a/7 — 1. 

As ^0 increases above /32, the solution changes its form. 
Instead of V having a minimum, the solution reaches a 
sonic point, where the equations do not determine the 
behaviour uniquely, so there can be a discontinuity [7]. 
In fact, only a subset of solutions are "regular" at a sonic 
point in the sense that the pressure gradient is finite and 
they can be extended beyond there. There are two possi- 
ble values of the pressure gradient, one associated with an 
isolated solution and the other with a 1-parameter fam- 
ily of solutions [41]. These regular transonic solutions do 
not contain black holes. Indeed, one can prove [3] that V 
cannot have a subsonic minimum and then rise through 
another sonic point to reach V = 1 again. Instead, they 
pass smoothly to the origin and represent either over- 
dense perturbations (for Bq < 0) or underdense pertur- 
bations (for ^0 > 0) which grow at the same rate as the 
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universe [7]. One can show that there is a continuum 
of underdense solutions and discrete bands of overdense 
ones [7, 41]. The solution with Bq = is exactly Fried- 
mann. These perturbed Friedmann solutions are of great 
interest in their own right [42] but they are not relevant 
to the present considerations. 

Since there is a one-parameter family of these solu- 
tions, we do not have to fine-tune the mass of the black 
hole to get self-similar growth as precisely as implied by 
Eq. (2.3). This is not surprising since the earlier analysis 
assumed that the accretion was driven by a local sound- 
wave rather than by an initial perturbation extending 
to infinity. However, the range of the parameter foi" 
which self-similar growth is possible is still narrow and 
this refiects the proximity of the Jean length and separate 
universe conditions. 

For comparison with the black hole solutions discussed 
later. Fig. 4 shows the Penrose diagrams for the acausal 
and universal black hole solutions. These correspond to 
Fig. 5 and Fig. 6 in Ref. [40], respectively. In these (and 
subsequent) Penrose diagrams a zig-zag line corresponds 
to a singularity, while a double-line corresponds to infin- 
ity. In the acausal case, it is worth emphasizing that the 
conformal structure is very different from that for a non- 
self-similar black hole in a Friedmann background. In 
particular, the black hole singularity is connected to the 
big bang singularity and necessarily naked for a while. In 
the universal case, the conformal diagram has two space- 
like singularities and there are no null infinities. This 
contrasts with the diagram for a recollapsing Fiedmann 
universe, which has two spacelike separated singularities. 




FIG. 4: The Penrose diagrams for an acausal black hole (top) 
and a universal black hole full (bottom) in an asymptotically 
quasi-Friedmann background. 



Later it was realized that these self-similar black hole 
solutions are not "properly" asymptotic to Friedmann 
after all [8]. The solid angle at infinity is no longer Air 
because the radial area - measured by the quantity S{z) 
- is perturbed out to arbitrarily large distances. More 
precisely, Eqs. (2.15) and (2.18) imply that the solid an- 



gle at infinity is 47re^^°, so there is a solid angle deficit 
in the black hole case {Bq < 0). These solutions may 
still be relevant to the real universe (since observations 
do not preclude such a deficit if it is not too large) but 
it would be more accurate to describe them as asymp- 
totically "quasi-Friedmann". Note that one tends to the 
separate-closed-universe condition as 5o ^ -co since 
the solid angle goes to zero. This is similar to the situ- 
ation with the Barriola-Vilenkin monopole, which is an 
approximate general relativistic solution for the triplet 
Higgs scalar field in the region far from the center [43]. 
The spacetime is asymptotically fiat in this monopole 
case but also contains a solid angle deficit. The prospects 
of detecting such a deficit by gravitational lensing obser- 
vations have been reviewed in Ref. [44]. 

The complete classification of single-fiuid shock-free 
positive-pressure spherically symmetric self-similar solu- 
tions by Carr and Coley [45] and Goliath et al. [46] has 
further clarified the situation. A key step in the Carr- 
Goley analysis is the derivation of all possible asymp- 
totic behaviors at large and small distances from the ori- 
gin. The behaviours at large spatial distances (usually 
corresponding to the limit z ^ oo) are of four kinds: 
(1) asymptotically quasi-Friedmann (1-parameter); (2) 
asymptotically quasi-static (2-parameter). There are also 
two families of solutions which exist only when 7 > 6/5: 
(3) asymptotically Minkowski at infinite z (1-parameter); 
and (4) asymptotically Minkowski at finite z but infinite 
physical distance (2-parameter). At small spatial dis- 
tances, the solutions are also of four kinds: they contain 
either (a) a black hole or (b) a naked singularity at finite 
z; or they can be connected to the origin at z = via a 
sonic point, in which case they are either (c) static or (d) 
a perturbation of the Friedmann solution there. 

The complete family of 1 < 7 < 2 solutions is found 
by combining the four kinds of large-distance behaviour 
and four kinds of small-distance behaviour. The way in 
which one connects the large-distance and small-distance 
solutions depends crucially on whether or not there is 
a sonic point. If the solutions remain supersonic ev- 
erywhere or subsonic everywhere, then the small- 2; be- 
haviour is uniquely determined by the large- 2; behaviour. 
The asymptotically quasi-Friedmann black hole solutions 
are everywhere supersonic and of type (la). However, if 
there is a sonic point, the solutions are more complicated 
because there can be a discontinuity. Indeed, we have 
seen that only a subset of solutions are regular at a sonic 
point. Although we focus here on the black hole asymp- 
totically quasi-Friedmann solutions, it should be stressed 
that these represent only a small subset of the total fam- 
ily of self-similar solutions. 

The study of critical phenomena in a Friedmann back- 
ground [47-51] may also be relevant. If one considers 
the evolution of primordial density perturbations with a 
variety of forms, whose amplitude is described by some 
parameter p, one finds that a black hole forms if p ex- 
ceeds some critical value and that the black hole mass 
scales as 

M oc{p-p^)^ (2.20) 
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for p ^ p^. The exponent /3 depends on the equation 
of state but is independent of the initial configuration of 
the cohapsing fluid. The critical solution itself is self- 
similar and has been identified explicitly for different 
values of 7 [52]. In terms of the Carr-Coley classifica- 
tion, it is either type (2d) (asymptotically quasi-static) 
for 1 < 7 < 1.28 or type (3d) or (4d) (asymptotically 
Minkowski) for 1.28 < 7 < 2. (Ah the critical solu- 
tions have a regular centre and need to be asymptotically 
Friedmann there). Although the critical solution is only 
attained in the limit for which the black hole mass goes to 
zero, the critical solution is only an intermediate attrac- 
tor, so there might still be a true self-similar attr actor in 
which the mass is non-zero. In any case, we have seen 
that there are at least asymptotically quasi- Friedmann 
black hole solutions with non-zero mass for this range of 
7- 

III. STIFF FLUID CASE 

For the limiting case in which 7 = 2, corresponding to 
a stiff fluid, the Zel'dovich-Novikov argument might con- 
ceivably apply because the sonic surface is a similarity 
horizon. Certainly this situation requires special consid- 
eration. Lin et al. at flrst claimed that there is a self- 
similar black hole solution in this case [9] . More precisely, 
they found a solution which has the exact Friedmann 
form outside the particle horizon (also a sonic point) but 
in which the velocity V inside this point has a minimum 
and then rises to 1 again, corresponding to a black hole 
event horizon. Such a solution resembles the ones shown 
in Fig. 3, except that it is subsonic between the event 
horizon and particle horizon. V therefore has a sub- 
sonic minimum, which is precisely the situation Carr and 
Hawking excluded in the 1 < 7 < 2 case. The pressure 
gradient appears to diverge at the event horizon but Lin 
et al. argued that this divergence can be removed by 
introducing an Eddington-Finkelstein-type coordinate. 

Later Bicknell and Henriksen [10] showed that this so- 
lution is invalid because the second point where V reaches 
1 as 2; decreases is not really a black hole event horizon. 
The reason for this is clarified in Ref. [15]. The metric 
induced on a const ant- 2; surface can be written as 

ds^ = -e^^(^)(l - V{zf)dt^ + r^S{zfdn^, (3.1) 

which is why V = 1 is usually a null surface. However, 
although Lin et al. pointed out that the pressure gradi- 
ent diverges where V = 1, they failed to appreciate that 
the density goes to zero there and this implies that e*^^^^ 
diverges. This means that the time component of the in- 
duced metric is non-zero, so the surface V = 1 is timelike 
rather than null and cannot be a black hole event horizon. 
However, it is worth stressing that Lin et al. correctly 
deduced the form of the solution up to the second point 
where V = 1. 

Bicknell and Henriksen did manage to construct a nu- 
merical self-similar stiff solution containing a PBH at- 
tached to an exact Friedmann exterior [10]. However, this 
solution has the feature that the stiff fiuid turns into null 



dust at the timelike surface where V = 1, with the fiuid 
becoming lightlike there. This situation might appear 
to be contrived but Reed and Henriksen [11] were able 
to find a solution of this kind by generalizing a previous 
solution of Hacyan [12] for a radiation-dominated model 
to arbitrary 7. In the Hacyan solution, the black hole is 
formed by infiowing radiation in a region described by the 
Vaidya metric. Reed and Henriksen showed that this so- 
lution is necessarily self-similar, with the black hole mass 
growing asymptotically as t, but it is generally unphysi- 
cal because it involves a decompression wave (which they 
claim violates the 2nd law of thermodynamics). However, 
they argued that it could be physical in the stiff case and 
that this might allow PBHs to form with masses up to 
4M0 if the equation of state is stiff up to the end of the 
hadron era. Although this seemed plausible at the time, 
the modern view does not favour a stiff hadron era, so 
this solution now also appears contrived. 

Recently, Harada et al. [15] have studied the stiff 
self-similar solutions with an exact Friedmann exterior 
in greater depth by treating the fiuid equations as a 
2-dimensional autonomous system and exploiting the 
equivalence (discussed below) between a stiff fiuid and 
a massless scalar field. They find a family of extensions 
within the inner V = 1 point which contain a massless 
scalar field. These might appear more natural than the 
Bicknell-Henriksen or Reed-Henrksen solutions but they 
still enter a negative mass region, which is presumably 
unphysical. 

Harada et al. also prove that there is no asymptotically 
Friedmann or asymptotically quasi-Friedmann solution 
with a sonic point containing a black hole in the stiff case. 
This is analogous to the situation in the non-stiff case, 
where solutions with sound-waves can only represent den- 
sity perturbations at the origin rather than black holes. 
The possibility of self-similar solutions is even more re- 
stricted in the stiff case because there is no range for the 
minimum of V between 1 and the sonic value. However, 
as discussed below, one might still expect the existence of 
universal self-similar black hole solutions (with Vmin > 1) 
in the stiff case, these being directly analogous to the su- 
personic solutions discussed in Section II. B. 



IV. SCALAR FIELD CASE 

Scalar fields are a key ingredient in modern cosmology. 
They play an important role during infiation and certain 
phase transitions; they feature in the preheating scenar- 
ios; and they are pervasive in all sorts of alternative the- 
ories (eg. string theory and scalar-tensor theory) which 
are likely to be relevant in the high curvature phase of the 
early universe. In such scenarios, scalar fields necessarily 
dominate the energy density of the universe at some stage 
of its evolution. One important physical feature in this 
case is that the scala r field can contain sound-waves, with 
the sound-speed {y^5p/5p 7^ ^Jpfp) being the speed of 
light in the short- wavelength limit [34]. Equation (2.12) 
then suggests that the accretion radius is \/2 times the 
Schwarzschild radius. As discussed later, this is different 
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from the case of a perfect fluid with negative pressure 
but constant 7, in which there are no sound-waves. 



A. Massless scalar field 

If there is no scalar potential, the scalar field is mass- 
less and one might expect the same results as in the stiff 
fluid analysis, since a scalar fleld is equivalent to a stiff 
fluid provided the gradient is everywhere timelike [13]. 
This is also consistent with numerical simulations of PBH 
growth in a Friedmann universe with a massless scalar 
fleld [14]. In fact, Harada et al. [15] have constructed 
an analytic proof that there is no self-similar PBH solu- 
tion in an exact or asymptotically Friedmann background 
if the universe is dominated by a massless scalar fleld. 
The proof exploits Brady's formulation of the self-similar 
equations as an autonomous system [56] and uses Bondi 
coordinates {v^r^O^cp) and the variables 



2M 

r 



(4.1) 



Here v is the ingoing radial null coordinate, M is the 
Misner- Sharp mass, =ln{r/\v\) and g = Qw/Ovr- 

The trajectories of the Harada et al. solutions in the 
(y, z) plane are illustrated in Fig. 5. The shaded regions 
are prohibited for a real scalar fleld. The gradient of the 
scalar fleld can change from timelike to spacelike on 7 = 
(thin solid line) or 22^(7 + /^) = 7 (thin dashed line) or 

= 0, where 



I'-- 



' 1 + 47r/^2 _ y- 
47r(l - 2z) 



(4.2) 



Here is a constant, which is taken to be l/\/127r. 
Two numerical solutions within the particle horizon are 
shown: the one above the Friedmann line is the posi- 
tive branch and goes directly to the negative mass region 
{y > 1), while the one below Friedmann reaches y = 3/4 
before going to the negative mass region. The Bicknell- 
Henriksen null-dust black hole solution connects to the 
lower numerical solution, but this is never realized by a 
scalar fleld. The black hole has an event horizon where 
the solution intersects z = 1/2 and an apparent horizon 
at y = 0, 2: = (X). 

Recent work by Bhattacharya et al [53] has found ex- 
plicit self-similar asymptotically Friedmann black hole 
solutions with Vmin > 1 in the massless scalar fleld 
case. These correspond to universal black hole solu- 
tions [40] and are a special subset of the Roberts self- 
similar solutions [54], the latter describing the inhomo- 
geneous spherically symmetric gravitational collapse of a 
massless scalar fleld which is minimally coupled to grav- 
ity (see also Ref. [55]). Such solutions should therefore 
also exist in the stiff case. 



B. Massive scalar field 

The situation is more complicated - but also more in- 
teresting - if there is a scalar potential (i.e. if the scalar 




FIG. 5: The massless scalar field solutions in the {y, z) plane. 
The expanding flat Friedmann solution is plotted as a thick 
solid line. Two numerical solutions within the particle hori- 
zon, whose significance is described in the text are shown by 
the thick dotted lines. The Bicknell-Henriksen null-dust black 
hole solution is indicated by a thin dashed-dotted line. 



field is no longer massless). In particular, this applies 
in the quintessence scenario, in which the scalar field 
rolls down a fiat potential [57]. In this context. Bean 
and Magueijo [19] have used a variant of the Zel'dovich- 
Novikov argument to claim that PBHs could grow up to 
~ IOOM0 through the accretion of quintessence before 
nucleosynthesis and this could be large enough to pro- 
vide the seeds for the supermassive black holes found in 
galactic nuclei. A generalization of this analysis by Cus- 
todio and Horvath [20], though not involving self-similar 
growth, has also claimed there could be appreciable ac- 
cretion of quintessence in some circumstances. 

We now examine these arguments in more detail. Fac- 
tors of c are included explicitly in this section for clarity 
and the dimensions of various quantities are justified in 
Appendix A. Bean and Magueijo [19] consider a scalar 
field with a potential of the form 



V{4>) = Voe- 



(4.3) 



where A is a dimensionless constant. This is the only form 
of the potential compatible with a self-similar solution, 
so one can immediately exclude the possibility of a black 
hole growing as fast as the universe with any other form. 
Note that (j) has units (? /\fG and V has units /{GL?) 
where L is a lengthscale. 

By generalizing an analysis of Jacobson [58] , which at- 
taches a Schwarzschild solution to a cosmological back- 
ground in which the field has the asymptotic value 0c (^) 5 
they claim that the energy fiux through the event horizon 

is (pc /c, where a dot denotes d/dt and t is cosmic time. 
This leads to a black hole accretion rate [59] 



dM 
~dt 



167rG^M^,^/c\ 



(4.4) 



They correctly note that only the kinetic energy of the 
scalar field is accreted. The scalar potential merely infiu- 
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ences the evolution of the background field, the potential 
and the cosmic scale factor: 



where is defined by 



c2(6 - A^) 



(4.5) 
(4.6) 



The expansion is accelerating providing A < y/2. Since a 
scalar field has 



(4.7) 



Eq. (4.4) can also be inferred from Eq. (2.9) of the 
Babichev et al. analysis if one assumes A = 4. From 
Eq. (4.4) and the expression for 0c given by Eq. (4.5), 
the accretion rate now becomes 



dM 



K 



SG 



(4.^ 



This equation can be integrated to give Eq. (2.3) but 
with the revised expression for K. 

As before, one has negligible growth for Mq <C K~^tQ, 
self-similar growth for Mq ~ K~^to and rapid divergence 
for Mq ^ K~^to. However, the self- similar mass is now 
of order the mass within the Hubble horizon (which is 
also the cosmological apparent horizon), since the parti- 
cle horizon does not exist if the acceleration extends back 
to the infinite past. More precisely, using the scale factor 
given by Eq. (4.5), the Hubble mass can be shown to be 



Mh 



IQttG 



3t 
2K 



(4.9) 



which is just a factor 3/2 larger than the mass corre- 
sponding to the self-similar solution. The self-similar 
result is the one which Bean and Maguiejo exploit. 
However, their analysis is exactly equivalent to that of 
Zel'dovich and Novikov [1], which Carr and Hawking 
proved was inapplicable. It is therefore important to 
know whether the non-existence of a black hole self- 
similar solution extends to this case. 

Custodio and Horvath [20] also consider black hole ac- 
cretion of a quintessence field with a scalar potential but, 
instead of Eq.(4.4), they use 



at 



(4.10) 



where the numerical factor differs because of relativistic 
beaming. They disagree with Bean and Magueijo's choice 
of the function 0c in Eq. (4.5) on the grounds that it 
neglects the local decrease in the background scalar field 
resulting from the accretion. For reasons which are not 
altogether clear, they therefore focus on a model in which 
the quintessence flux into the black hole is constant. This 
requires that the potential have a specific form, which is 



different from (4.3) and therefore incompatible with self- 
similarity. This leads to the mass evolution 



M: 



Mo 



K'Mo 
to 



it - to) 



(4.11) 



where is a constant related to the (fixed) flux. This is 
equivalent to the Zel'dovich-Novikov formula (2.3) pro- 
viding one replaces K hy K't. In particular, the mass 
diverges at a time 



to 1 



K'Mo 



(4.12) 



and they attribute this unphysical feature to the fact that 
the constant-fiux assumption must fail. They therefore 
consider alternative models in which the flux decreases as 
a power of time. In general, they find that the increase 
in the black hole mass is small unless its initial value is 
finely tuned. 

Although Custodio and Horvath conclude that the 
Bean-Magueijo argument is wrong, we would claim that 
they have not identified the more significant problem - 
that it neglects the background cosmological expansion. 
Indeed, since the cosmological expansion is neglected in 
both these analyses, they could both be flawed for the 
same reason as the original Zel'dovich-Novikov one. 

To examine the relativistic situation, we seek a spher- 
ically symmetric self-similar solution which represents 
a black hole in a quintessence-dominated cosmological 
background. In the quintessence scenario, for which the 
similarity assumption requires the potential to have the 
form given by Eq. (4.3), the analysis in Ref. [15] shows 
that there is no self-similar black hole solution in an ex- 
act or asymptotically Friedmann or asymptotically quasi- 
Friedmann background if the universe is decelerating and 
no such solution in an exact Friedmann background if it is 
accelerating. However, this does not prove non-existence 
for the case in which the universe is accelerating and 
asymptotically Friedmann, which is what would be re- 
quired to disprove the Bean-Magueijo proposal. 

Kyo et al. [21] have shown that, with an accelerat- 
ing scalar potential, there is a one-parameter family of 
self-similar solutions which are properly asymptotically 
Friedmann, in the sense that there is no solid angle 
deficit. In such solutions, the perturbation falls off very 
rapidly at spatial infinity. The question of whether such 
solutions can have a black hole event horizon remains 
open. However, as discussed in the next section, there 
is a self-similar black hole solution in a closely-related 
scenario. 



V. DARK ENERGY FLUIDS 

The observed acceleration of the universe means that 
violation of the strong energy condition is required. The 
definition of the various energy conditions is given in Ap- 
pendix B; for a perfect fluid with ^> = (7 — l)pc^ and 
positive p, the violation of the strong one corresponds to 
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7 < 2/3 (including negative values). We describe such 
a fluid as "dark energy" . However, it should be stressed 
that this term is used in various ways by different authors 
and indeed is sometimes used interchangeably with the 
term "quintessence" , so some clarification of terminology 
is required. In general, a scalar field corresponds to a 
stiff fluid (7 = 2) when the potential is negligible and 
to a cosmological constant (7 = 0) when the potential 
dominates, so it corresponds to a fluid with < 7 < 2 at 
intermediate times and to "dark energy" (in our sense of 
the term) for any period when < 7 < 2/3. Quintessence 
corresponds to the special case in which the scalar field 
rolls down a very flat potential, so that 7 may change 
but is always in the latter range [57]. We have seen that 
there are sound-waves in the scalar field case but there 
are no so und-w aves for a perfect fluid with < 7 < 1 
because \ppj~p = ^7 ~ 1 c is imaginary. 

The question of whether a PBH can grow in a self- 
similar manner if the universe is dominated by dark en- 
ergy can be addressed using the analysis of Babichev et 
al. [65]. Although the parameter A given by Eq. (2.10) 
is not well-defined for < 7 < 1, Babychev et al. ar- 
gue that it can be taken to be 4 throughout this range. 
In this case, Eqs. (2.2) and (2.11) can still be applied, 
which suggests that self-similar growth is again possible. 
However, note that ^ 00 as 7 ^ 0, so the self-similar 
black hole has zero mass in this limit. 

As in the positive-pressure case, one should be wary 
of the self-similar claim since the above analysis neglects 
the cosmic expansion. The exact relativistic analysis has 
been carried out in our two recent studies of self-similar 
solutions with dark energy [22, 23]. These use a com- 
bination of numerical and analytical methods to clas- 
sify all spherically symmetric self-similar solutions with 
< 7 < 2/3 which are asymptotically Friedmann at large 
distances. This extends the Carr-Coley classification of 
self-similar solutions with 1 < 7 < 2 to the < 7 < 2/3 
case. In fact, Ref. [60] already includes an asymptotic 
analysis of such solutions, although there are some errors 
in that work. 

There are several differences from the 1 < 7 < 2 sit- 
uation. First, the fact that there are no sound-waves 
changes - and indeed simplifies - the analysis, since there 
can be no discontinuities and solutions are analytic ev- 
erywhere. Second, the limiting values of z for large 
and small spatial distances are reversed in the Fried- 
mann solution: large spatial distances now correspond 
to 2: ^ and small ones correspond to 2: ^ cxd. Fur- 
thermore, we will find that some solutions necessarily 
span both positive and negative z values. There are 
also several differences as regards the solutions them- 
selves. For example, one loses the exact static solution 
ifO < 7 < 2/3 because one needs positive pressure 
to balance gravity. This also means that there are no 
asymptotically static ones, although there are still what 
are termed asymptotically quasi-static solutions (a point 
missed in Ref. [60]). On the other hand, one gains some 
solutions, since the Kantowski-Sachs and asymptotically 
Kant owski- Sachs solutions [61] now become physical. Al- 
though one loses the solutions for 6/5 < 7 < 2 which 



are asymptotically Minkowski at large distances [45] , one 
gains solutions which are asymptotic to a negative-mass 
Schwarzchild singularity at small distances. 

The < 7 < 2/3 case is particularly important in 
the cosmological context because there then exists a one- 
parameter family of self-similar solutions which are prop- 
erly asymptotic to the fiat Friedmann model at large dis- 
tances. This has also been pointed out in Ref. [62]. As 
we have seen, this does not apply in the positive pres- 
sure case because the solutions are then only asymptoti- 
cally quasi- Friedmann. However, these solutions are still 
acausal, in the sense that there is no sonic point, so they 
are not like the solutions which Carr and Hawking were 
originally seeking in the positive pressure context. 

In order to investigate these solutions, we have used 
the asymptotic analysis of Ref. [22] and integrated the 
field equations numerically in Ref. [23] . We find that the 
solutions which are asymptotically Friedmann at large 
distances {z 0) are described by a single parameter Aq, 
which is analogous to the parameter Bq in the 1 < 7 < 2 
case and may be interpreted as the density perturbation 
at spatial infinity. However, the expressions for A and 
B in Eqs (2.18) and (2.19) no longer apply and must be 
replaced by 

A ^ A^z^^-^^l^, B ^ -(Ao/67)^^'"^^/^ (5.1) 

in the < 7 < 2/3 regime. The form of the solutions 
changes as one varies the parameter and one has three 
possible behaviours [23] . The first class of solutions (with 
positive) is confined to the z > domain and con- 
tains a naked singularity; the Friedmann solution itself 
has = 0. The second class (with A^ negative but 
not too small) involves an extension into the 2: < re- 
gion and contains a black hole; this is analogous to the 
Carr-Hawking solution in the positive-pressure case. The 
third class (with A^ even more negative) also involves an 
extension into the z < region and corresponds to a 
wormhole solution; this has no analogy in the positive- 
pressure case. 

The form of V in these solutions is shown as a function 
of 1/2; in Fig. 6 for the 7 = 1/3 case and we discuss them 
below in more detail. Note that V = 1 now corresponds 
to a cosmological event horizon rather than a particle 
horizon. Although we did not investigate all values of 
7 in the range to 2/3, we would expect the results 
to be qualitatively similar. We note that the system of 
equations was also numerically investigated in Ref. [62] 
but no black hole or wormhole solutions were reported 
there because the analytic extension of the solutions to 
the 2: < domain was not considered. 

A. Black hole solutions 

Solutions with ai < Aq < for some critical 7- 
dependent value ai (i.e. which are not too overdense) are 
asymptotically quasi-Kant owski- Sachs as 2: ^ oc. The 
exact Kantowski-Sachs metric has the form 

^ _ (2 - 37) (2 - 7) ^^2 ^^4(l-7)/7rf^2 ^^2^^2 (5 2) 

T 
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2 I 1 '.'!'!! ' r~| 0.02 

I 1 i i : 0.01 

i i i i : 

• ! ! ! / -0.01 

; I j j .-; -0.02 

1.5 - i i i i i / - -0.03 




• Naked singularity 
Friedmann 

• Black hole 

• Wormhole 



FIG. 6: The form of V in the negative pressure solutions 
with 7 = 1/3, showing the transition from the black hole to 
wormhole to naked singularity solutions as one varies the pa- 
rameter Aq which describes their asymptotic behaviour. The 
ordinate is taken to be —l/z, which is large and negative in 
the asymptotic Friedmann region and small near the origin if 
the solution reaches there. 



with a suitable radial coordinate and this means that 
the physical distance tends to a finite limit (like a shell) 
rather than zero as z ^ +oo. However, these solutions 
can be extended analytically into the region with nega- 
tive z and, as z increases from — oo, they then approach a 
positive-mass singularity at some negative value z^ . Since 
the Kantowski-Sachs solution has a curvature singularity 
at t = 0, the analytic extension in this case must be in- 
terpreted as an extension from the positive r to negative 
r region. However, the areal radius and mass are still 
positive in the extended region. This solution describes 
a black hole in an asymptotically Friedmann background 
and is therefore analogous to the Carr-Hawking solution 
in the positive pressure case. In particular, there are 
two similarity horizons, corresponding to a cosmological 
event horizon dit z = Zi and a black hole event horizon 
at z = 2^2, and a black hole singularity at z = z*. How- 
ever, the initial big bang singularity is null rather than 
spacelike in this case. The Penrose diagram for this so- 
lution is shown in Fig. 7 and is reproduced from Fig. 12 
of Ref. [23]. 



z=+0 




FIG. 7: The Penrose diagram for a solution with negative 
pressure which contains a black hole in an asymptotically 
Friedmann universe. 

For 7 = 1/3, we have found numerically that the ratio 
of the size of the black hole event horizon to the Hubble 



length goes from and 0.70 as Aq decreases from to 
ai. This means that the size of a self-similar cosmological 
black hole has an upper limit, so only sufficiently small 
black holes can grow as fast as the universe, but may be 
arbitrarily small. Although a black hole may not grow 
as fast as the universe when it first forms (if it is as 
large as the background Hubble length), the ratio of its 
radius relative to the Hubble length will decrease as time 
proceeds. Eventually it will fall below 0.70, after which 
the black hole starts to grow self-similarly. 

Since there is a one-parameter family of these solu- 
tions, we do not have to fine-tune the mass of the black 
hole in the way implied by Eq. (2.3) to get self-similar 
growth. This also applied in the positive pressure case 
but the parameter range is now even more extended. Also 
the positive-pressure solutions may be less plausible since 
they are only asymptotically quasi- Friedmann. 

B. Wormhole and white hole solutions 

Solutions with as < Aq < ai (where as is another 7- 
dependent critical value) are also asymptotically quasi- 
K ant owski- Sachs as z ^ 00. These solutions can again 
be analytically continued into the negative z region but 
they now reach z = 0~ (corresponding to infinite phys- 
ical distance) instead of a singularity. They are asymp- 
totically Friedmann at z = 0+ but asymptotically quasi- 
Friedmann at 2; = 0~, except for a special value a2 be- 
tween ai and 0^3 which allows the solution to be asymp- 
totically Friedmann at both ends. Since the physical ra- 
dius never becomes zero except at the big bang singular- 
ity, these solutions correspond to cosmological wormholes 
and there are two cosmological event horizons at zi and 
^2. The Penrose diagram for these solutions is shown in 
Fig. 8 and is reproduced from Fig. 15 of Ref. [23]. This 
only differs from Fig. 7 in that there is no singularity at 
the top left. 
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FIG. 8: The Penrose diagram for solutions with negative 
pressure in which a wormhole connects an asymptotically 
Friedmann solution on the right to an asymptotically quasi- 
Friedmann universe on the left. 

The transition from the black hole to wormhole solu- 
tions has a simple physical explanation. As the param- 
eter Aq becomes more negative, the ratio of the size of 
the black hole apparent horizon to the size of the cosmo- 
logical apparent horizon increases. The transition occurs 
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when this ratio reaches unity and thereafter both hori- 
zons disappear, every region being trapped. This is in 
contrast to the 1 < 7 < 2 case, where the two apparent 
horizons never merge and one tends to a separate closed 
universe as the black hole size increases [5]. 

The solutions with Aq < as are asymptotically quasi- 
static rather than asymptotically quasi- K ant owski- Sachs 
as 2; ^ 00 but they can still be extended into the neg- 
ative z region. These solutions have a central curvature 
singularity at some negative value z*, where the mass is 
positive, and correspond to the formation of a wormhole 
from a white hole in a Friedmann background at t = 0. 
The Penrose diagram for this solution is shown in Fig. 9 
and is reproduced from Fig. 15 of Ref. [23]. 



z=+0 




FIG. 9: The Penrose diagram for a solution with negative 
pressure which describes the formation of a wormhole from a 
white hole in an asymptotically Friedmann spacetime. 

In order to understand these wormhole solutions, some 
clarification of what is meant by a wormhole is required. 
This is discussed in detail in another paper [63] but we 
summarize the main points here. A wormhole is an ob- 
ject in general relativity which connects two or more 
asymptotic regions. The most famous example of a static 
wormhole is the Morris- Thorne solution [66] and it is 
well known that this requires a violation of the null en- 
ergy condition [67, 68]. The definition of this is given 
in Appendix B but it corresponds to /i + p > for a 
perfect fiuid. Dynamical wormholes are not so well un- 
derstood but their study was pioneered by Hochberg and 
Visser [69] and Hay ward [70], who defined a wormhole 
throat in a dynamical spacetime as some kind of trap- 
ping horizon [71]. In their definitions, a wormhole throat 
is a two-dimensional surface of non-vanishing minimal 
area on a null hypersurface and the null energy condition 
must again be violated there. However, there is no past 
null infinity in the context of our cosmological wormhole 
solutions because there is an initial singularity. In fact, 
the Hochberg- Visser and Hayward definitions are inappli- 
cable because the spacetime is trapped everywhere and 
there is no trapping horizon. This demonstrates that 
these definitions miss the important class of cosmologi- 
cal wormholes which are asymptotically Friedmann and 
start with a big bang. 

In order to remedy this problem, we define a worm- 
hole throat as a two-dimensional surface of non- vanishing 



minimal area on a spacelike hypersurface [63]. The one- 
parameter family of spherically symmetric self-similar 
wormhole solutions in an accelerating Friedmann back- 
ground discussed above satisfy this definition. They are 
asymptotically Friedmann at one infinity and they have 
another infinity, which may also be asymptotically Fried- 
mann for a special value of the parameter. Interestingly, 
the dominant energy condition is satisfied everywhere. 

In Ref. [63] we construct two analytic examples of self- 
similar cosmological wormholes. One corresponds to the 
numerical solution obtained above, but it contains a sin- 
gular hypersurface which violates the null energy condi- 
tion. The other is a smooth model involving a combi- 
nation of a perfect fiuid and a ghost scalar field (i.e. a 
massless scalar field with a negative kinetic term) but 
the total matter content still satisfies the dominant en- 
ergy condition. 

VI. EXOTIC CASES 

In this section, we consider values of 7 which are out- 
side the "conventional" range from to 2 but might nev- 
ertheless arise in more exotic scenarios. These cases have 
not been studied in sufficient detail to come to definite 
conclusions. In particular, the spherically symmetric self- 
similar solutions have not been obtained explicitly. Nev- 
ertheless, it is worthwhile summarizing what is known 
about these situations. 



A. Negative pressure fluids satisfying SEC 

A fiuid with 2/3 < 7 < 1 has negative pressure but 
it is not sufficiently negative to violate the strong en- 
ergy condition. Self-similar solutions in this context have 
been studied in Ref. [60] , which shows that their possible 
large and small distance behaviours include asymptoti- 
cally quasi-Friedmann and asymptotically Kasner (cor- 
responding to a black hole singularity), respectively. In- 
deed, the equations are formally similar to those in the 
1 < 7 < 2 case, except for the absence of sound-waves 
and that should not be relevant if one is seekng an ana- 
logue of the 1 < 7 < 2 solutions which are supersonic 
everywhere. One would therefore expect the black hole 
solutions to exist, although this has not yet been demon- 
strated numerically. 

B. Phantom fluids 

A fiuid with 7 < is described as a "phantom" . More 
generally, this corresponds to any fiuid with p < —p and 
this possibility has been explored in the context of fiu- 
ids with equation of state of the form p = a{p — po) 
in Ref. [64]. Here we confine attention to the situation 
with a = 7 — 1 and po = because of the similarly 
assumption. There are still asymptotically Friedmann 
self-similar solutions in this case, so the previous analy- 
sis can in principle be extended to cover this. However, 
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since the scale factor is a{t) = there is a "big 

rip" singularity [73, 74] at t = 0, spacelike infinity is at 
z = +00 and the domain of definition is t < 0. Note 
that there is no future null infinity in these solutions, so 
there is no black hole event horizon. However, there is 
still a black hole apparent horizon and one may take this 
to define the mass of the black hole. Equation (2.9) then 
has the consequence that the PBH mass should decrease 
as a result of phantom accretion and it would need to 
decrease in proportion to the time left to the big rip in a 
self-similar solution. 

Babichev et at. [65] have analysed black hole accre- 
tion of phantom energy and find that the mass evolves 
according to 

M= ^, ^ Mo 1^1, (6.1) 



where Mi is the initial mass and Mq = {'^^c^ /2GA) with 
A taken to be 4. The first expression uses the Babichev 
et al. time coordinate with the time of the big rip cor- 
responding io t' = r; the second expression applies in 
the limit t' ^ r and uses our time coordinate t. Note 
that Eq. (6.1) is very similar to Eq. (2.3) with K be- 
ing given by Eq. (2.11). Although the mass evolution 
is clearly self-similar close to the big rip, the derivation 
of Eq. (6.1) again fails to account for the background ex- 
pansion, since the black hole is not embedded in a cosmo- 
logical background. Nevertheless, self-similar black hole 
solutions might in principle exist in this situation. The 
fact that black hole mass is shrinking rather than increas- 
ing is of obvious astrophysical interest. There have also 
been studies of black hole accretion in a ghost conden- 
sate [72, 75]. However, this term should not be confused 
with the word "ghost" which is sometimes used for the 
phantom fiuid itself. A ghost condensate is more like dust 
(7=1) than a phantom fiuid. 

C. Tachyonic fluids 

A fiuid with 7 > 2 might be described as "tachyonic" , 
in the sense that the sound-speed exceeds c. It formally 
contradicts the dominant energy condition but not the 
other ones. In fact, most of the equations in the 1 < 
7 < 2 case be formally extended into this domain. In 
particular, all the curves in Fig. 2 can be extended to 
arbitrarily large 7. Babichev et al. [80] have considered 
black hole accretion in this case but the problem does 
not yet seem to have been studied in the cosmological 
context. 



VII. DISCUSSION AND SUMMARY 

In this review, we have discussed the possible existence 
of self-similar solutions containing a black hole or worm- 
hole in an asymptotically Friedmann background whose 
density is dominated by a perfect fluid with p = (7— l)pc^ 



or a scalar field. A simple Bondi-type analysis predicts 
an accretion rate of the form given by Eq. (2.2) in all 
cases, though with a different value of the constant 
which appears to permit self-similar growth. Ultimately, 
this is because the Friedmann equation implies that the 
density scales as p oc in a fiat Friedmann background, 
which is precisely the condition for self-similarity. 

However, the simple analysis is suspect because it ne- 
glects the cosmic expansion, so this has motivated a 
more careful relativistic analysis which allows for the 
expansion. In the positive-pressure case, the Bicknell- 
Henriksen solution (in which a stiff fiuid turns into null 
dust) seems to be the only known self-similar black hole 
solution with a sound-wave which is exactly Friedmann at 
large distances and even this is rather contrived. There is 
also the Hacyan radiation-dominated solution (in which 
the region containing the black hole is described by a 
Vaidya solution) but this is also artificial because all the 
photons have to become radially directed within some 
point and it also violates the 2nd law of thermodynam- 
ics. 

Apart from these examples, it appears that there is 
no self-similar solution containing a black hole in either 
an exact or asymptotically Friedmann background for 
any value of 7 in the range 1 < 7 < 2. There are 
only asymptotically quasi-Friedmann solutions, includ- 
ing "universal" black holes without a black hole event 
horizon or cosmological particle horizon. However, there 
are self-similar asymptotically Friedmann solutions for 
< 7 < 2/3, which suggests that PBHs can grow as fast 
as the universe in the presence of dark energy (at least for 
a limited period). This conclusion may also apply for a 
quintessence field, although this has not been rigorously 
proved. 

The difference between the positive and negative pres- 
sure solutions is important in two respects. First, while 
the negative-pressure ones are physically well- motivated 
in the infiationary scenario, because one might expect the 
associated density perturbations to extend to "infinity" , 
the positive-pressure ones are theoretically unmotivated 
and may also be observationally excluded for some pa- 
rameter because they have a solid angle deficit at large 
distances. Second, self-similar black holes only exist if 
their size as a fraction of the cosmological horizon is not 
too small in the positive-pressure case but not too large 
in the negative-pressure case. This means that less fine- 
tuning is required in the latter case. It is interesting that 
there is no accretion in the limit 7^0, which is consis- 
tent with the Schwarzschild-de Sitter solution. 

The existence of these self-similar black hole solu- 
tions suggests that black holes may increase their mass 
by a considerable factor during any dark-energy or 
quintessence dominated era, regardless of whether or not 
they are "primordial" [19] . There are two contexts in 
which this effect may be important: (1) in the period 
immediately after any PBHs formed at the end of infla- 
tion; (2) in the recent period when whatever dominates 
the density of the universe causes it to accelerate. Bean 
and Magueijo [19] focussed on the flrst situation but the 
second one may also be interesting. For although the 
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accretion factor may not be very large up to now, the 
black hole mass will continue to grow like cosmic time so 
as long as the dark energy dominates the cosmological 
density. In the simplest models, this applies indefinitely, 
so the black hole can grow arbitrarily large. We are in- 
vestigating the astrophysical implications of this result 
- and especially its implications for the Bean-Magueijo 
claim that PBHs can grow large enough to provide the 
supermassive black holes in galactic nuclei - in a separate 
paper. 

Finally, we note that our analysis restricts the situa- 
tions in which the similarity hypothesis [76] applies. This 
hypothesis claims that there are certain circumstances 
in which spherically symmetric solutions evolve to self- 
similar form. The present work shows that there are 
at least some situations in which this does not happen. 
(Even if there were a self-similar solution, one would 
still need to show that it was stable in order for it to 
be an attractor.) On the other hand, the similarity hy- 
pothesis does appear to hold in spherical gravitational 
collapse when the pressure is positive but very small 
(0 < 7 — 1 <C 1) [77, 78]. It may also hold in the negative- 
pressure situation. 
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Appendix A: Physical dimensions 



We take the action to be 



/ 



dtd^xy/—g 



16nG 



n 



(Al) 

where the coordinates in 4-dimensional Lorentzian space- 
time are (x^, x*) = (ct, x*) with i = 1, 2, 3. We emphasize 
that dt d^x must not be written as d'^x if we include fac- 



tors of c explicitly. The Einstein equations are 



(A2) 



V = - i^9^.v<)^A'' j - 9f.uV. (A3) 

The dimensions of the action and gravitational constant 
are 

[S]=ML^T-^, (A4) 

[G] = M-^L^T-^, (A5) 

so we obtain 

[cVG] =MLT-^, (A6) 

[dtd^x X cVG] =ML^T-\ (A7) 

and hence 

[T^,]=ML-^T-\ (A8) 

The relativistic factor in the black hole metric is 

[GM/c^] = L. (A9) 

For a perfect fluid, the energy momentum tensor is 

T^u = {fi+p)u^'u^ -^pg^,, (AlO) 

where = dx^/d{cT) is the dimensionless 4- velocity 
of the fluid element, with r being the affine time. The 
mass density p, energy density /i and pressure p have 
dimensions 



so 



[p] =ML-\ 
[11] =ML-^T-^, 
[p] =ML-^T-^, 

[Gn/c'] = [Gp/c^] = [Gp] = 
For a scalar field, we have 

[V] =[c^/{GL^)] = ML-^T-^. 



(All) 
(A12) 
(A13) 

(A14) 

(A15) 
(A16) 



Hence the exponential potential can be written in the 
dimensionally correct form 



(A17) 



where A is a dimensionless constant. The equation of 
motion for 6 is 



dV 



(A18) 



In the cosmological context, a scalar field is equivalent to 
a perfect fluid with 



2\dx^ ) 



1 



(A19) 
(A20) 
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For the flat Ftiedmann metric 



3' 



where is the fiat line element, and 

rr.t 



27rGA4t2 • 



(A21) 



(A22) 



This is the energy density for a Friedmann universe filled 
by a scalar field with an exponential potential. The last 
two equations also apply for a perfect fluid with p = 
(7 — l)/i if one replaces A with \/^. 



Appendix B: Energy conditions 

We summarize the energy conditions for a matter field 
with energy-momentum tensor given in the diagonal form 
as T^j, = diag(— /i,Pr,Pt,Pt) [81, 82]. The physical inter- 
pretations of /i, Pj. and pt are the energy density, radial 
pressure and tangential pressure, respectively. 

The inequality T^yW^W^ > for any timelike vec- 
tor is called the weak energy condition (WEC). 
This guarantees that a timelike observer measures a non- 
negative energy density and implies > 0, Pr + M ^ 



and Pt + M ^ 0- The null energy condition (NEC) just re- 
places with a future-directed null vector in the above 
definition, which means that one can drop the condition 
/i > 0. 

The condition that T^^Wi, is a future-directed non- 
spacelike vector for every future-directed non-spacelike 
vector is called the dominant energy condition 

(DEC). This means that the mass-energy can never 
flow faster than the speed of light and implies /i > 0, 
— M ^ < M —fi < Pi < fi. 

The inequality (T^^ - g^^T^ p/2)W^W > for any 
timelike vector is called the strong energy condi- 
tion (SEC). If there is no cosmological constant, this as- 
sures the timelike convergence condition R^yW^W^ > 
in general relativity, which means that the gravitational 
force is essentially attractive. The SEC implies Pr-h/i > 0, 
Pt + M > and pr + 2pr + /i > 0. 

For = Pt = Pi as assumed throughout this paper, 
the WEC imphes /i > and p + /i > 0, the DEC implies 
/i > and — /i < p < /i, the SEC implies p -\- > and 
3p + /i > 0, and the NEC implies /i + p > 0. The third 
and fourth conditions must be violated for inflation and 
the dynamical wormholes, respectively. Note that DEC 
implies WEC, WEC implies NEC, and SEC implies NEC, 
but SEC does not imply WEC. 
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